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Abstract. Combinatorial expressions are presented to the solutions to initial value problems 
of the discrete and ultradiscrete Toda molecules. For the discrete Toda molecule, a subtraction- 
free expression of the solution is derived in terms of non-intersecting paths, for which two 
results in combinatorics, Flajolet's interpretation of continued fractions and Gessel-Viennot's 
lemma on determinants, are applied. By ultradiscretizing the subtraction-free expression, the 
solution to the ultradiscrete Toda molecule is obtained. It is finally shown that the initial value 
problem of the ultradiscrete Toda molecule is exactly solved in terms of shortest paths on a 
specific graph. 



PACS numbers: 02.30.1k, 04.20.Ex, 05.45.Yv 

1. Introduction 

The Toda molecule HI is a semi-infinite version of the Toda lattice 12 ■ As is the Toda lattice, 
the Toda molecule is known as a typical example of integrable systems which has, despite its 
nonlinearity, an exact solution in terms of Hankel determinants. The discrete Toda molecule 
is a discrete analogue of the Toda molecule which is derived in (3) by using the bilinear 
formalism. The discrete Toda molecule is a discrete integrable system which possesses a 
Hankel determinant solution analogous to the Toda molecule. 

The discrete Toda molecule is also derived by using the Lax formalism 0|5), in which a 
connection with orthogonal polynomials is exploited to deduce the time evolution equations 
of the discrete Toda molecule 

In + e n ~ In + e „+V ^ l& > 
OH) (r) 0) (lb) 

ef = forfeZ, « e N„. (lc) 

In numerical algorithms, the equations (fl3 are known as recurrence equations of the qd 
algorithm, which is used for computing Pade approximants of analytic functions (see, e.g., 
||6l), and for computing eigenvalues of tridiagonal matrices (see, e.g., Q). In the study 
of pure combinatorics, Viennot |8] applied the qd algorithm to a combinatorial problem 
of enumerating configurations of non-intersecting paths. From the viewpoint of dynamical 
systems, Viennot's combinatorial result is observed as solving an initial value problem of the 
discrete Toda molecule (Q3 from particular initial value. 

In recent progress of investigating integrable systems, much attention has been given to 
ultradiscrete integrable systems, especially since the discovery of a direct connection between 
discrete integrable systems and soliton cellular automata (9). A general method to derive an 
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ultradiscrete integrable system from a discrete integrable system is called ultradiscretization. 
For the discrete Toda molecule (|T), ultradiscretization is performed as follows IfTUl : Introduce 
new dependent variables Q„, E„ by q„ - exp(-(2„ f) /e), e„ = exp(-£ , ^' ) /s) with a parameter 
s > 0, and take the limit s — > 0. The equations ([TJ then tend to the ultradiscrete Toda molecule 



(t l) =™{£(g-Z4? 1 \l%X (2a) 

I k=0 k=0 ) 

Ci" = Qti ~ Qt +l) + E nli - for ? 6 Z, « e No. (2b) 
It is shown in IfTUl that the ultradiscrete Toda molecule (O describes the dynamics of a box- 
ball system. 

In this paper, we examine the initial value problems of the discrete and ultradiscrete 
Toda molecules, ([T]i and (|2j, for which the initial value is given at t = 0. Obviously, one 
can exactly solve the problems in the following sense: At any time t e No and any site 
n e No, the exact value of each dependent variable can be calculated from the initial value 
in finitely many arithmetic and minimizing operations. However, it is still nontrivial how 
to formulate the solutions since the equations are nonlinear. The aim of this paper is to 
derive an exact expression of the solutions to the initial value problems purely in terms of 
the initial value. In order to formulate the solutions, we will utilize combinatorial objects, 
non-intersecting paths and shortest paths on a graph, in view of the combinatorial results 
on paths: Flajolet's interpretation of continued fractions ifTTl and Gessel-Viennot's lemma 
lfT2l[T3l on determinants. 

This paper is organized as follows. In section|2] we review a determinant solution to the 
discrete Toda molecule, based on which, in section [3] we combinatorially formulate an exact 
expression of the solution to the initial value problem of the discrete Toda molecule in terms 
of non-intersecting paths. In section |4] we derive the solution to the initial value problem of 
the ultradiscrete Toda molecule by ultradiscretizing the solution to the discrete Toda molecule 
obtained in section [3] Further combinatorial observations lead us to a simpler expression of 
the solution in terms of shortest paths on a specific graph. Sectionals devoted to concluding 
remarks. 



2. Determinant solution to the discrete Toda molecule 



In section [2] we give a brief review on a determinant solution to the discrete Toda lattice 
together with bilinear equations associated with the discrete and ultradiscrete Toda molecules. 
See, e.g., lfT4l [Kfl for detailed explanations. Based on the determinant solution, in the 
subsequent sections, we will examine initial value problems of the discrete and ultradiscrete 
Toda molecules. 

We introduce a tau function r„ of the discrete Toda molecule ([]} by the variable 
transformation 

_(f+i)_(0 Jf+VJt) 
a® = " +1 " e m = " n+1 77 6 No (3) 

'n+1 n+1 n+1 

for which we assume the boundary condition that r? = 1 . We then obtain from ([TJ a bilinear 
equation of the discrete Toda molecule, 

'n+1 'n+1 'n ' „+2 ^ 'n+1 'n+1' « t 
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To the bilinear equation (0), we have an exact solution in the Hankel determinant of size n 
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where /„ is an arbitrary function subject to the linear dispersion relation 

(The determinant of size zero is assume to be unity conventionally.) We can verify the solution 
© by means of Sylvester's determinant identity. Substituting the determinant (Bat to (01, we 
obtain an exact solution to the discrete Toda molecule (0. 

We can derive a bilinear equation of the ultradiscrete Toda molecule (J2J by 
ultradiscretizing the procedure for the discrete Toda molecule: Introducing a tau function 



(5b) 



7f by 
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with r*" = 0, we then obtain from (f2]i a bilinear equation of the ultradiscrete Toda molecule, 



1 n+\ n+\ 



min{r ( ( / +1) + r"; 9 1) ,2r ( '? l ) 

1 « n+2 ' n+l J 



n e N . (7) 

The equations ©, (O for the ultradiscrete Toda molecule are obtained by 
ultradiscretizing the corresponding (]3), (0} for the discrete Toda molecule: Assume that 
qf = expi-Q^/e), ef = exp(-£ , J/ ) /e), = ex.p(-T®/e) with s > and take the limit 
e — > 0. Then, the equations ©, (H]i tend to ©, (O, respectively. The limiting procedure in 
ultradiscretization replaces the operations x, /, + with +, -, min, respectively, for we have 

- e log(e~ x/£ x e- Y/s ) =X + Y, (8a) 

-s log(e- x/e /e- Y/E ) =X-Y, (8b) 

-elog(e^ /e + e~ Y/E ) min {X, Y] as s -> 0. (8c) 

However, the counterpart of subtraction, -, does not exist since the limit of s log(e~ x / £ - e^ E ) 
as s — > is undetermined. Commonly, we cannot ultradiscretize equations containing 
subtractions. Especially, we cannot directly ultradiscretize the determinant solution (O 
because we may encounter subtractions in expanding the determinant. In section[3] we derive 
a subtraction-free expression of the determinant solution (0 to which ultradiscretization is 
directly applicable. 



3. Initial value problem of the discrete Toda molecule 

As an initial value problem of the discrete Toda molecule, we consider the following: 
For the discrete Toda molecule (0, let us write the initial value at t — 

9 f = fl2 „, e^=a 2n+ u «€N . (9) 

Then, for every t e No, find the exact value of 'q„ \ e v n e No, uniquely determined from ((T) 
in terms of the initial value a„. 

To the qd algorithm for Pade approximants whose recurrence equations are given by (JTJ, 
a combinatorial interpretation was given by Viennot 0, in which a combinatorial expression 
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Figure 1. A positive grounded path P of steps UUDDUUUDDUDD. The label of each step, 
a„ or 1, is shown below the step. The path P weighs w(P) = 0^0^02- 



of the determinant (Bat is formulated in terms of non-intersecting paths. The fundamental 
idea used in section[3]comes from Viennot's approach to the qd algorithm. 

The discrete Toda molecule ([TJ is linearized in the following sense: The nonlinear system 
(Q~|) for q„ , e„ reduces into the linear system ( |5bl for fjf' through the dependent variable 
transformations © and ( |5at via the tau function r„ . We can thus evaluate the time evolution 
of the discrete Toda molecule by the dispersion relation (l5bl whose initial value problem is 
exactly solved by 

j®=j™, t,n€ N . (10) 
The initial value problem therefore amounts to the following two subproblems: 

(i) Find the initial value fif^ of /„ (,) at t = in terms of a„ from (0, (01 and (|5}. 

(ii) Find the value of the determinant Oat to evaluate the tau function t„ for each t e No. 

In order to solve the subproblem dj}, we utilize Flajolet's combinatorial interpretation of 
continued fractions IfTTI : Let us consider a path P in Z 2 consisting of up steps U = (1,1) and 
down steps D = (1, -1). We say that P is positive if P never goes beneath the x-axis, y = 0. 
(A positive path can touch the x-axis.) We say that P is grounded if both the two ends, the 
initial and terminal points, of P are on the x-axis. Figure Q] shows an example of a positive 
grounded path P. We label each step in P by a n if the step is an up step ascending from the 
line y — n and by unity if a down step. We then define the weight w(P) of P by the product of 
the labels of all the steps in P. For example, the path P in figure Q] weighs w(P) = a^a 3 ^. We 
conventionally assume the weight of empty paths of no steps to be unity. We refer by D(P) to 
the number of down steps in P. 

Lemma 1 (Flajolet UTTl ). It holds that 
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«<^" = -aoz <"> 
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1 

where the (formal) sum in the left-hand side is taken over all the positive grounded paths P 
whose initial point is fixed at (0, 0). 

The subproblem (01 asks us to solve the system of equations 

- A " +lA " « 2 „ +1 = ^%^ «eN„, (12a) 



rm 

A„ = det(^»)-i , A,', = det(/^ +1 )-l . (12b) 



for /„ (0) , n e No, where A„ and A', denote the determinants of size n 
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In the theory of Pade approximants (see, e.g., ||6]), it is well-known that the system ( IT2b is 
equivalent to the equation between a formal power series and an S -fraction 

»=° 1 - aiz 
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with the normalization that = 1 . 

This observation on Pade approximants leads us to the solution to the subproblem (01: 
Owing to lemmaQ] with the normalization that = 1, 

/f = (14) 

p 

where the sum in the right-hand side is taken over all the positive grounded paths P whose 
two ends are fixed at (0, 0) and (2n, 0). For example, the first few of /„ (0) are 

(15a) 
(15b) 

a\, (15c) 
'qAi + floOj + a(,a\a2- (15d) 
As noted in IfTTI . we can write the combinatorial formula (fT4l in the form 

k x + \k 2 + \ + l 

fn } = XI Yj Y " ' Y ak < aklClh ' ' ' ak » ■ ( 16) 

*i=Ok=o* 3 =o k„=a 

The initial value /„ <0) of /„ (,) is thus found as a polynomial in a^ homogeneous of degree n. The 
number of monomials in /„ (0) , which is equal to the number of positive grounded paths from 
(0,0) to (2«,0), is counted by the Catalan number C„ = ^( 2 n "). For details on the Catalan 
numbers, refer the On-Line Encyclopedia of Integer Sequences 1151 for Sequence A000108. 

Due to dTOb and ( fPfl i, the (7, £)-entry of the determinant (J5aj>, f® k = f^- +k , is shown 
to be equal to the sum of the weight w(P) of all the positive grounded paths P going from 
(0, 0) to (2t + 2j + 2k, 0), or equivalently, going from (-2 j, 0) to {It + 2k, 0). We can thereby 
successfully apply Gessel-Viennot's lemma on determinants lfl2l[T3l to solve the subproblem 
(juj: For t, n € No, let P(t, n) denote the collection of n-sets P = {Pq, . . ., P n -\ } of positive 
grounded paths Pj satisfying the following conditions: 

(a) The two ends of Pi are fixed at (-2 j, 0) and (2t + 2 j, 0). 

(b) The « paths Pq, . . ., P n -\ are non-intersecting, Namely, every two distinct paths Pj and 
Pk, j + k, never intersect at any points. 

Then, with the normalization that - 1, Gessel-Viennot's lemma yields that 

if = w W where w(P) = w(Po)---w(P„- 1 ). (17) 

PeP(t,n) 

As shown in figure |2] we can draw each n-set P = {Pq, . . ., P n -\ } € P(t, n) as a diagram of n 
positive grounded paths which are non-intersecting. 

The value of r„ is found as a polynomial in a^ homogeneous of degree n(2t + n — l)/2. 
The number of monomials in t„ ', which is equal to the cardinality #P(t, n), is exactly evaluated 
in II, 

#Pit , n)= FT 2 JL1211. (18) 

l<J<k<t J 
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Figure 2. The n-sets P = {Pq, . . ., P n -i | e P(t, n) of positive grounded paths, where t = 4, 
n = 3. Each n-set P can be drawn in a diagram of n non-intersecting positive grounded paths 
Pj such that Pj goes from (-2 j, 0) to (2t + 2j, 0). 



We have solved the subproblems (01 and The solution to the initial value problem of 
the discrete Toda molecule is given as follows: 

Theorem 2. The solution to the initial value problem of the discrete Toda molecule (fJ]) is 
given by with the tau function 

rf= w ( p )> t,n€N . (19) 

PeP(i.n) 

The expression (IT~9T > of the tau function r„ is subtraction-free, namely, contains no 
subtractions. That is why t„ is positive for every t, n e No if and only if the initial value a„ is 
positive for every n 6 No- The subtraction-free expression (TT9) of t„ is ultradiscretizable to 
obtain an exact solution to the ultradiscrete Toda molecule. 

4. Initial value problem of the ultradiscrete Toda molecule 

As an initial value problem of the ultradiscrete Toda molecule (0, we consider the totally 
analogous problem to the discrete Toda molecule solved in section[3j 

For the ultradiscrete Toda molecule (0, let us write the initial value at t — 

= A 2n , fiW = A2n+1> „ e No _ ( 20) 

r/zen, /or every f 6 No, /«^/ f/ie exacf vaZwe o/ Q„ , v n e No, uniquely determined from 
dU) in fermi of the initial value A n . 

The solution to this initial value problem can be obtained by ultradiscretizing the 
corresponding solution to the discrete Toda molecule stated in theorem |2] Ultradiscretizing 
theorem |2] we obtain the following statement: 

Let P be a positive grounded path. We label each step in P by A„ if the step is an up 
step ascending from the line y — n and by zero if a down step. We define the weight W(P) 
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Figure 3. The hooks H in a positive grounded path P. The path P contains three up hooks, 
highlighted in dotted lines, and two down hooks, in dashed lines. 

by the sum of the labels of all the steps in P. The solution to the initial value problem of the 
ultradiscrete Toda molecule is then given by ((6) with the tau function 

T®= min W(P), W(P) = W(P ) + --- + W(P„-i). (21) 

PeP(t,n) 

The support P(t, n) of the minimum is much the same as the sum in ( fl9] l. 

The rest of section|4]is devoted to simplifying the expression d2"TT i of the tau function T„ 
by combinatorial observation. 

4.1. Solution in tabular paths 

Let P be a positive grounded path. We refer to two consecutive up-down steps and down- 
up steps, UD and DU, in P by a peak and a valley, respectively. We say that P is tabular 
provided that, for some k e No, all the peaks and the valleys in P reside in the strip of height 
one bordered by the two horizontal lines y = k and y = k + 1. 

For t, n e No, we define a subset P{t,n) c P(t,n) as the collection of n-sets P = 
{Pq, . . . , P n -\ } e P{t, n) in which every Pj is tabular. For example, all the n-sets P e P{t, n) in 
figure |2 except the upper right one, belong to P(t, n) since each path in P is tabular. 

Lemma 3. There exists P e P(t, n) which takes the minimal weight: W(P) — mmp^, ^ W(P). 

Proof. We will prove the lemma by examining particular subpaths which we call hooks. Let 
P be a positive grounded path. We call a subpath H of P an up hook (resp. a down hook) 
provided that H is of at least four steps and that both the first and the last steps of H are 
down steps (resp. up steps) and all the middle steps are up steps (resp. down steps). That is, 
each up hook (resp. down hook) is of the form DU k D (resp. UD k U) for some integer k > 2. 
For example, see figure [3] Obviously, a positive grounded path P is tabular if and only if P 
contains no hooks. 

Let us define two maps (p and if/ which deform a positive grounded path P as follows: 
ip{P) (resp. if/(P)) denotes the positive grounded path obtained from P by replacing each up 
hook in P, say DU k D, with U k l DUD (resp. with DUDU k l ), where <p(P) = \f/(P) = P if P 
contains no up hooks. See figure [4] which shows the deformed paths <p(P) and if/iP) obtained 
from the positive grounded path P in figure [3] 

We can observe the following on the maps ip and \p: Let P and P' be positive grounded 
paths. 

(a) The ip and if/ never increase the number of up hooks in P as well as the number of down 
hooks in P. 

(b) If P contains no up hooks to the right side of the line x — j then tp(P) so to x — j - 2. 
Similarly, if P contains no up hooks to the left side of x — j then if/{P) so to x = j + 2. 

(c) The <p and if/ never move the location of the two ends of P. 
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Figure 4. The deformation by the maps ip and tjj. The paths <f(P) and if/(P) (depicted in solid 
lines) are obtained from the path P (in dotted lines) in figure[5]by deforming each up hook in 
P as indicated by an arrow. 



(d) If P and P' are non-intersecting then tp(P) and <p(P') are so. That is also the case with 
if/(P) and t//(P'). 

(e) The following 'mean' formula holds: 

W((p(P)) + WW)) = 2W(P). (22) 

Now, let us prove lemma [3] Let P = {Pq, . . . ,P n -i] e P(t,n). We can assume that P 
takes the minimal weight, W(P) = min/» e f(, „) W(P'), without any loss of generality. For 
each k e No, let if k (P) = {tp k (Po), • ■ ■ , <p k (P n -i)} denote the «-sets of positive grounded paths 
obtained from P by applying the map (p iteratively k times to each path Pj e P. By induction 
with respect to k e No, we can show the following: 

(i) Due to the observation ©, every path in (p k (P) contains no up hooks to the right side 
of the line x = 2t - 2k - 1 . That is because P contains no up hooks to the right side of 
x = 2t- 1. 

(ii) Due to the observations (JcJ and (0, tp k (P) € P(t, n). 

(iii) Due to the observation fle), W(<p k (P)) = W(P). Indeed, if W((p k (P)) > W(<p k - l (P)), the 
formula (l22l would lead W{i[f o tp k ~ l (P)) < W(tp k ~ l (P)), that contradicts the minimality 
of W(P). 

As a consequence of (i), (ii), (iii), we can deduce that (p'~ l (P) 6 P(t, n) contains no up hooks 
and has the minimal weight W(tp'~ l {P)) = W(P). 

In a similar way, we can show the existence of P e Pit, n) containing no down hooks 
and having the minimal weight W(P) = W(P). From the discussion in the last paragraph, the 
tp'~ l (P) contains no up hooks and has the minimal weight W((p (P)) = W(P). Further, due 
to the observation ©, the (p'~ l (P) contains no down hooks. Therefore, the ip'~ l (P) having the 
minimal weight belongs to the set P{t, n). That completes the proof. □ 

As a consequence of lemma [3] we can also solve the initial value problem of the 
ultradiscrete Toda molecule in the following way: 
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Figure 5. The graph G and a path Q e Q(t, n), where t = 4, n = 3. The graph G splits into 
two disjoint subgraphs, one of which is drawn in solid lines and the other in dashed lines. The 
path Q, drawn in thick lines, goes between t) and (; + 2n, 1). 



Theorem 4. The solution to the initial value problem of the ultradiscrete Toda molecule d2]) 
is given by (|6|) with the tau function 

T® = min W(P), (,neN . (23) 

PeP(r,n) 

The expression d23l of the tau function T„ is supported by the set P(t, n) much smaller 
than P(t, n) in OTT i. In fact, the cardinarity of P(t, n) is equal to the binomial number 

(t + n — l\ i — r n + j 
• 1 II <24) 

which is much smaller than #P(t, n) given by (IT8T l. In this sense, the expression ( 1231 gives a 
simpler expression of the tau function 7*„ than (fJTJ. 



4.2. Solution in shortest paths on a graph 

In section [4~2l based on theorem|4] we derive another combinatorial expression of the solution 
in terms of shortest paths on a graph. 

Let G denote the (directed acyclic) graph in Ny consisting of the vertices at the points 
(j, k) e Nq, j > k, connected by the two types of (directed) edges: east edges E ^ = (2, 0) and 
south edges S j,k = (0, -1). Here the subscripts j,k indicate that the initial points of Ej^ and 
S j,k are at the vertex (j, k). As shown in figure|5] since the east edges Ej,k have length two, the 
graph G splits into two disjoint subgraphs. 



Combinatorial solutions to the discrete and ultradiscrete Toda molecules 



10 



We define the weight function W over the edges in G by 

j-k-l 

W(E j>k ) = J] At + kAj- k , W(S U ) = 0, (25) 
e=o 

where A„ are the initial value for the ultradiscrete Toda molecule. We think of the weight 
W{e) of an edge e as the length of e. For each path Q on G, we then think of the weight W(Q) 
as the length of Q, which is equal to the sum of the weight W(e) of all the edges e passed by 
Q. (The length W(Q) may be negative due to the arbitrariness in A n .) 

For t e No, t > 1, and n e No, let Q(t,n) denote the collection of paths on G between 
the two vertices (f, f) and (f + 2n, 1). We then have a one-to-one correspondence between 
P = {Po, ■ • ■ ,Pn-i] e P(t,n) and Q e Q(t,n): The tabular positive grounded path Pj has its 
peaks and valleys in the strip bordered by the lines y = k - 2j and y = k - 2j + 1 if and only if 
the path Q on G passes through the east edge E t +2j,t-k- For example, in the «-sets P e P(t, ri) 
of positive grounded paths in figure[2] the lower left one belongs to P(t, n) and is in one-to-one 
correspondence with the path Q in figure [5] Actually, the weight function W on G is defined 
in (l25T l so that W(P) = W(Q) for every pair of P and Q in one-to-one correspondence. We 
thereby have the identity 

min W(Q) = min W(P). (26) 

QeQ(t,n) PeP(t,n) 

For /, n e No, let us define Q{t, n) to be the collection of paths on G between the two 
vertices (f, f) and (f + 2n, 0). We can then show that the identity d26l i still holds even if we 
replace the support set Q(t, «) of the left-hand minimum with Q(t, «), 

min W(Q) = min W(P), (27) 

QeQd.n) PeP(t,n) 

for W(Ejfi) = W{Ej y \) for every j e No, j > 1. In addition, in that case, the identity ( |27T ) also 
takes place for t = 0. 

Finally, combining theorem[4]and the identity ( TZTI i. we obtain the following result: 

Theorem 5. The solution to the initial value problem of the ultradiscrete Toda molecule (0 
is given by ([6} with the tau function 

r, ( ,° = min W(Q), f,neN , (28) 

QeQ(t,n) 

where the weight function W on the graph G is defined by ( I25l >. 

The right-hand side of (l28l denotes the length of the shortest paths on the graph G 
between the two vertices (t, f) and (t + 2n, 0). It should be noted that Nakata lfl6l constructed 
a similar combinatorial expression in terms of shortest paths (called minimum weight flows 
in ifTBI ) of a particular solution to the ultradiscrete Toda molecule on the finite lattice, 
n = 0, 

5. Concluding remarks 

In this paper, we have investigated the discrete and ultradiscrete Toda molecules from a 
combinatorial viewpoint. To the tau function which solves an initial value problem of 
the discrete Toda molecule, we have given a combinatorial expression in terms of non- 
intersecting paths. Especially, in order to read the tau function in combinatorial words, 
we utilized Flajolet's path interpretation of continued fractions and Gessel-Viennot's lemma 
on determinants and non-intersecting paths. Due to the combinatorial expression, we have 
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succeeded to derive a subtraction-free expression of the tau function which is given in a 
Hankel determinant. 

For an initial value problem of the ultradiscrete Toda molecule, we first obtained an exact 
solution by ultradiscretizing the tau function of the discrete Toda molecule. We next rewrote 
the tau function to derive a simpler expression in terms of shortest paths. As a result, we have 
shown that the tau function which solves the initial value problem of the ultradiscrete Toda 
molecule can be evaluated as the length of shortest paths on a specific graph in which the 
length of edges is determined by the initial value. 

In this paper, we deduced combinatorial expressions of the tau functions with the help of 
a determinant solution to the discrete Toda molecule and the technique of ultradiscretization. 
For the tau functions ( fT9] i and (l28l l given in terms of combinatorial objects, however, it is 
expected to make combinatorial (or bijective) proofs to directly verify that the tau functions 
satisfy the bilinear equations © and Q. For that purpose, the technique of alternating walks 
IfTTll for Schur symmetric functions would be useful. 

The combinatorial idea used in this paper should be applicable to other discrete integrable 
systems associated with continued fractions, such as the R\ and Z?n chains with R\- and in- 
fractions UHlQjD, the FST chain with the Thiele-type continued fractions ||201 . and the matrix 
qd algorithm with the matrix S-fractions lED . Those applications will be discussed in future 
works. 
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